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INTRODUCTION 


The  US  Navy  currently  uses  or  is  planning  to  use  a  number  of  electrooptica!  (EOl 
systems.  l  or  example,  forward-looking  infrared  (El.IR)  thermal  imaging  systems  are 
currently  installed  on  many  aircraft  as  an  integral  part  of  their  weapon  systems.  Infrared 
search  and  track  (IRS  I  )  systems  are  currently  in  a  research  and  development  stage.  Heat- 
seeking  or  laser-guided  weapons  have  been  used  for  a  long  time  and  new  systems  are  under 
development.  To  optimize  the  performance  of  these  FO  systems  in  tactical  situations,  a  way 
to  assess  their  performance  is  needed.  For  surveillance  systems,  this  means  a  way  to  predict 
the  sy  stem's  ability  to  detect  and  correctly  identify  a  target  of  interest;  for  weapon  systems, 
a  way  to  predict  the  system's  ability  to  track  and  destroy  the  target. 

An  assessment  of  surveillance  and  tracking  equipment  naturally  requires  a  knowl¬ 
edge  of  the  limitations  and  capabilities  of  the  operator  and  sensor  hardware,  but  it  also 
requires  a  know  ledge  of  the  thermal  characteristics  of  the  target  and  its  natural  back¬ 
ground.  and  the  transmissivity  between  the  target  and  sensor.  These  latter  parameters  are 
required  its  inputs  to  computer  algorithms  w  hich  determine  the  ranges  at  which  the  equiv¬ 
alent  blackbody  temperature  difference  between  the  target  at  zero  range  and  its  back¬ 
ground  is  degraded  by  the  atmospheric  transmittance  to  the  minimum  detectable  or  resolv¬ 
able  temperatures  of  the  sensor. 

Both  target  and  background  infrared  signatures  are  controlled  by  the  intervening 
atmosphere.  For  a  surface  ship,  heating  from  the  sun  and  internal  sources,  and  cooling  by 
the  ambient  winds  and  ship's  motion  (as  well  as  solar  insulation  effects)  cause  the  ship's 
temperature  to  vary.  The  apparent  temperature  difference  between  the  ship  and  its  back¬ 
ground.  as  viewed  by  the  sensor,  will  then  depend  upon  the  infrared  transmittance  and 
emission  along  the  propagation  path.  To  be  of  practical  use  in  tactical  situations,  airborne 
FI  1R  systems  must  be  able  to  detect  and  identify  surface  targets  at  ranges  of  several  tens  of 
kilometers.  This  requirement  limits  the  depression  viewing  angles  to  within  a  few  degrees  of 
the  horizontal.  As  an  example,  the  range  to  a  target  as  viewed  from  an  altitude  of  1  km 
with  a  I -degree  depression  angle  is  80  km.  Depressing  the  viewing  angle  to  3  degrees  places 
the  target  only  20  km  away.  For  larger  viewing  angles,  the  infrared  radiance  of  the  :;ca 
surface  is  the  primary  contributor  to  the  background  scene.  Closer  to  the  horizon,  however, 
the  sky  radiance  reflections  and  emissions  by  the  intervening  atmosphere  must  be  taken 
into  account.  While  a  smooth  sea  is  a  poor  emitter  of  infrared  radiation  at  the  small  view¬ 
ing  angles,  emissions  from  the  indiv  idual  wave  facets  of  a  wind-ruffled  sea  contribute  to  the 
background  scene. 

I  he  development  of  suitable  target  and  background  models  to  be  used  in  computer 
algorithms  of  sy  stem  performance  presents  a  formidable  task  and  is  needed  by  both  users 
and  designers  ol  FO  systems.  In  1978.  the  Under  Secretary  of  Defense  for  Research  and 
Engineering  formulated  goals  to  be  followed  by  the  three  services  to  address  the  primary 
factors  (e.g..  aerosol  extinction,  gaseous  absorption,  refraction,  and  turbulence)  affecting 
i he  transmission  and  emission  characteristics  of  the  atmosphere  which  determine  sy  stem 
performance.  Basically.  these  goals  (as  updated  in  1982)  are 

•  Accvratclv  model  propagation  environments  of  naturally  occurring  and  man¬ 
made  aerosols  and  gases,  including  the  effects  of  turbulence  and  multiple  scat¬ 
tering.  Relate  the  particulate  constituents  to  measurable  and  predictable  meteo¬ 
rological  parameters. 

•  Develop  atmospheric  sensor  systems  lor  model  validation,  test,  and  operational 
support.  Special  emphasis  should  be  given  to  studies  of  the  applicability  of  sens¬ 
ing  by  light  detecting  and  ranging  (lidar). 
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•  Relate  atmospheric  effects  to  the  performance  of  military  systems,  in  opera¬ 
tional  as  well  as  in  research  and  development  activities,  to  include  the  develop¬ 
ment  of  tactieal  decision  aids.  This  includes  the  consideration  of  transmission, 
background,  and  target  factors  which  are  based  on  climatologies  or  real-time 
en\ ironmental  inputs. 

Since  their  inception,  considerable  progress  has  been  made  toward  meeting  these 
goals.  The  effects  of  gaseous  absolutions  on  radiance  and  transmittance  have  been 
adequately  addressed  in  the  LOWTRAN  7  code  developed  at  the  Air  force  Geophysical 
I  aboratory  (Knei/ys  et  ai..  1988).  Since  multiple  scattering  by  aerosols  has  recently  been 
included  in  I  OW7  RAN  7  and  turbulence  has  been  addressed  in  the  Army’s  Electro-Optical 
Systems  Atmospheric  Effects  1  ibrary  (EOSEAE)  (Miller  and  Ricklin,  1987),  they  wall  not 
be  covered  here.  In  the  following  sections,  examples  will  be  presented  of  advances  made 
toward  meeting  each  of  the  goals,  with  primary  emphasis  on  the  measurement  and  model¬ 
ing  of  absorption  and  scattering  by  aerosols  and  their  effects  on  radiance  and  transmit¬ 
tance.  w  hich  determine  system  performance.  It  is  the  purpose  of  this  document  to  summa- 
rve  past  accomplishments,  discuss  present  deficiencies,  and  define  future  directions  for  the 
assessment  of  Navy  EO  system  performance. 


MARINE  AEROSOL  MODELING 


Selectable  aerosol  si/e  distribution  models  are  now  available  in  LOWTRAN  7  for 
calculating  their  scattering  and  absorption  properties.  Of  primary  interest  here  is  the  Navy 
Maritime  Aerosol  Model  (NAM)  (Gathman,  1983).  This  model  was  developed  by  using  an 
extensive  set  of  measurements  made  in  the  Eastern  Pacific  (San  Nicolas  Island)  during 
1978.  1979.  and  1980.  These  measurements  were  of  broadband  (visible  through  the  far 
infrared)  atmospheric  transmission  and  meteorological  parameters,  aerosol  size  distribu¬ 
tions,  and  radon  concentrations.  This  was  a  cooperative  effort  between  Navy  Laboratories 
(Naval  Ocean  Systems  Center.  Naval  Research  Laboratory,  Naval  Surface  Warfare  Center, 
and  Naval  Weapons  Center),  the  Pacific  Missile  Test  Center,  and  the  Naval  Postgraduate 
School.  Shipboard  measurements  in  the  Atlantic  Ocean  (USNS  Hayes ,  1977;  USNS  Lynch, 
1983)  of  aerosol  si/e  distributions  and  meteorological  parameters  were  also  included  in  the 
model’s  development. 

Key  to  the  development  of  NAM  was  the  recognition  of  a  trimodal  structure  in 
aerosol  si/e  distributions  plotted  from  the  aforementioned  field  experiments.  This  trimodal 
distribution  was  interpreted  as  a  superposition  of  three  aerosol  distributions;  the  smallest 
due  to  continental  aerosols,  an  intermediate  si/e  dependent  on  the  average  wind  speed,  and 
the  largest  governed  by  the  instantaneous  wind.  Accordingly,  a  si/e  distribution  model  (at 
radius  r)  is  the  sum  of  three  log-normal  distributions  given  by 
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Component  -1 ,  represents  the  contribution  by  continental  aerosols.  AM  is  an  air 
mass  parameter  that  is  allowed  to  range  between  integer  values  ol  1  tor  open  ocean  and  10 
tor  coastal  areas  and  is  given  by 

AM-  IX'KRn  4)  +  I  (5) 

where  Rn  is  the  measured  atmospheric  radon  content  expressed  in  pC'i  m'  In  the  absence 
of  radon  measurements,  the  ait  mass  factor  can  be  related  to  the  elapsed  time.  7~ (days),  tor 
the  air  mass  to  reach  the  point  of  observation. 

.1.1/  =  /.Y7'[9  expt  7  4;]+  I  (6) 

Components  .-t ,  and  .1 ,  represent  equilibrium  sea  spray  particles  generated  by  the  surface 
wind  speed  averaged  over  24  hours  (  l  .  in  m  s)  and  the  current  surface  wind  speed  (  C  in 
m  s).  respectively.  In  Eu.  I.  rr  the  modal  radius  for  each  component  referenced  to  a  rela¬ 
tive  humidity  of  80C  (r ,  =  0.03  pm.  r:  =  0.24  pm.  and  rx-  2.0  pm),  is  allowed  to  grow  with 
relativ e  humidity  (  RH)  according  to  the  Fitzgerald  ( 1975)  formula: 

F-  [(2  RH  100)  6<  I  RH  100)] 1  '  (7) 

The  contribution  to  the  total  extinction  or  absorption  by  each  component  can  then  be 
w  ritten  as 


o,.M  (A),  =  I SF )  |  C,  /  Qt.JK.r.m )  exp 
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w here  C,  =  (0.00 1  n!  f)A  r  The  factor /  1  in  the  expression  for  C,  ensures  a  constant  total 
number  of  particles  us  the  relative  humiditv  increases.  )  is  the  cross-section  for 

either  the  extinction  or  absorption  normalized  to  the  geometrical  cross-section  of  the 
spherical  particle,  and  m  is  the  complex  refractive  index,  which  is  allowed  to  change  from 
that  of  dry  sea  salt  as  the  particle  deliquesces  with  increasing  humidity.  LOWTRAN  6 
provides  precalculated  values  in  tabular  form  of  the  parameter  ol,a(Kj)  Ct  at  discrete  wave¬ 
lengths  tor  four  relative  humidif  es  (50f7.  85f7.  90ri .  and  99r().  front  which  the  average 
extinction  for  a  specific  wavelength  band  and  relative  humidity  can  be  readily  determined 
by  interpolation.  When  an  observed  surface  v  isibilitv  ( l  IS,, )  is  av  ailable  as  an  input  to  the 
model,  the  amplitudes  of  the  three  components  are  adjusted  by  a  scaling  factor  (SF)  so  that 
the  calculated  aerosol  extinction  coefficient,  a, .  at  a  wavelength  of  0.55  pm.  is  the  same  as 
the  observed  extinction,  n  determined  from  the  relationship 


I  IS. 


3.912 


a  .  +  a. 


(9) 


w  here  n r  is  the  Rav  leigh  contribution  to  extinction  at  0.55  pm. 

7  his  model  has  undergone  extensive  evaluation  of  the  accuracy  with  w  hich  it  can 
predict  infrared  radiance  and  transmission.  It  should  be  noted  that  the  current  wind  speed 
component.  A ,.  is  different  from  the  value  published  in  I.OWTR  AN  6  ( Knei/ys  et  al.. 
1983).  which  was  Ax  =  0.01527(7  2.2).  This  modification  was  found  to  be  necessary  in 

order  to  match  published  measurements  of  infrared  (8  12  mm)  sky  radiances  ( Hughes. 
1987)  In  that  work,  it  was  determined  that  for  low  wind  speed,  the  measured  and  calcu¬ 
lated  radiances  agreed  within  2f7  at  the  optical  horizon.  However,  during  moderate  wind 


speed,  the  current  wind  speed  component  of  the  aerosol  model  had  to  be  lowered  by  factors 
near  0.05  for  the  calculated  and  measured  radiances  to  agree.  When  the  model  was  used, 
the  discrepancies  between  measured  and  calculated  infrared  radiances  were  found  to  be 
insensitive  to  the  24-hour  averaged  wind  component  (Fig.  1).  Measurements  by  de  Leeuw 
( 1986)  of  large  particle  si/e  distributions  (r  >  5  microns)  provided  another  method  of 
evaluating  the  current  wind  speed  component.  In  that  work,  size  distributions  were 
measured  with  an  impactor  at  different  altitudes  (0.2  to  1  I  meters)  above  the  sea  surface.  In 
order  tor  the  size  distributions  calculated  with  the  model  using  measured  surface  meteoro¬ 
logical  parameters  to  agree  with  the  measured  values,  the  current  wind  speed  component 
had  to  be  lowered  by  a  factor  near  0.07  (Fig.  2).  Within  the  measurement  accuracies  of 
both  the  radiance  and  particle  size  comparisons,  the  reduction  factors  in  the  A ,  component 
were  in  reasonable  agreement. 

F he  accuracy  to  which  the  model  can  predict  extinction  coefficients  has  also  been 
tested  against  the  transmissometer  and  meteorological  measurements  at  San  Nicolas  Island 
(Gathman  and  L  iters.  1983).  Good  correlations  were  obtained  between  calculated  and 
measured  extinctions  for  near  and  midinfrared  wavelengths.  However,  for  the  far  infrared, 
the  calculated  extinctions  were  20rJ  to  400  greater  than  those  measured  by  the  transmis¬ 
someter.  While  these  correlations  may  have  been  influenced  by  the  uncorrected  Ay  compo¬ 
nent.  they  were  shown  to  be  sensitive  to  the  selection  of  the  air  mass  factor  and  to  w  hether 
the  visibility  was  used  as  an  input. 

The  model  was  dev  eloped  to  be  as  representative  as  possible  of  different  atmos¬ 
pheric  conditions:  however,  it  cannot  be  expected  to  exactly  reproduce  the  optical  proper¬ 
ties  in  a  given  location  at  any  specific  time.  Neither  of  the  techniques  for  selecting  the  air 
mass  factor  (radon  count  or  air  trajectory  analysis)  is  currently  available  for  shipboard  use. 
Also  the  visibilities  inferred  from  point-scattering  devices  onboard  ship  are  most  apt  to  be 
contaminated  by  ship  effluences.  A  method  is  needed  for  selecting  the  input  parameters  so 
that  the  model  best  represents  a  particular  situation.  A  remote-sensing  technique  has 
recently  been  developed  ( Hughes  and  Jensen.  1988)  whereby  unique  values  of  the  air  mass 
factor  and  visibility  can  be  selected  for  different  meteorological  conditions.  These  values 
are  inferred  from  l.OWTRAN  calculations  which  agree  with  both  the  surface  measure¬ 
ments  of  8-12-/zm  horizon  radiances  and  visible  atmospheric  optical  depths  determined 
from  satellite-detected  upwelling  solar  radiances. 

For  this  study,  a  Piper  Navajo  aircraft,  equipped  with  Rosemount  temperature  and 
pressure  probes  and  an  EG&G  dewpoint  sensor,  made  vertical  spirals  over  the  ocean  to 
obtain  the  profile  of  temperature,  relative  humidity,  and  pressure  which  are  required  inputs 
to  the  LOWTR AN  6  computer  code  for  calculating  the  sky  radiances.  (A  Barnes  PRT-5 
radiation  thermometer  was  also  onboard  the  aircraft  to  measure  sea  surface  temperatures 
from  low-level  constant-altitude  flights.)  At  the  time  the  meteorological  parameters  were 
obtained  on  29  September  and  25  November  1987.  measurements  of  infrared  (812  pm) 
horizon  radiances  were  also  made  with  a  calibrated  thermal  imaging  system  (AG A 
I  HERMOVISION.  model  780)  having  a  2.95-degree  field-of-view  lens.  The  thermal  video 
processor  system  (THER  MOTE  KNTX)  available  with  the  AGA  system  displays  the  ther¬ 
mal  scene  on  a  computer  screen  in  a  format  consisting  of  128  pixel  lines  (0.023  deg  pixel 
line).  The  response  of  the  system  is  determined  by  placing  a  blackbodv  of  known  tempera¬ 
ture  (  n().  I°C  for  temperatures  <50°C)  in  front  of  the  lens  aperture.  The  digitized  v  ideo 
signal  transfer  function  of  the  system  then  allows  the  blackbodv  temperature  to  be  repro¬ 
duced  to  within  ±0.2'' C. 
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vc,  CURRENT  WIND  SPEED  (m/s) 

Figure  1.  Sensitivity  of  the  LOWTRAN  6  NAM  to  surface  wind  speed 
parameters. 


RADIUS  Uim) 

Figure  2.  Comparisons  of 
measured  particle  size  distri 
butions  with  distributions 
calculated  with  the  original 
and  adjusted  LOWTRAN  6 


I  or  these  measurements,  the  scanner  was  located  at  an  elevation  of  33  meters  on 
the  Point  1  tuna  peninsula  in  San  Diego  and  was  directed  due  west  over  the  ocean  so  that 
appt o.ximatelv  halt  the  field-of-v  iew  was  above  and  half  below  the  horizon.  The  data 
processing  software  of  the  AG  \  system  also  allows  the  the  effective  biackbody  temperature 
til  each  pixel  in  the  scene  to  be  determined.  On  each  of  the  days,  the  pixel  corresponding  to 
the  maximum  temperature  in  the  thermogram  (20.4°C  on  29  September  and  I6.5°C  on  25 
November)  was  taken  to  coincide  with  the  infrared  horizon.  From  the  current  and  24-hour 
averaged  wind  speeds  (  P  =  7.5  m  s  and  (  =  4.4  m  s  for  29  September,  and  P.  =  2.5  m  s 
mid  I  -  3.9  in  s  lor  25  November)  measured  on  shore  at  the  AGA  site,  and  the  vertical 
proliles  o!  meteorological  parameters.  1  OW  1  RAN  6  calculations  were  made  to  agree  with 
the  maximum  pixel  radiance  in  the  scene  hv  using  nonunique  combinations  of  air  mass 
lactois  and  \  isibilities. 

In  these  calculations,  the  meteorological  profiles  were  divided  into  33  lasers,  as 
allowed  hv  1  ()\V  I  R  AN  6  I  he  lower  layers  of  the  profiles  are  also  div  ided  into  sublayers 
containing  the  same  amount  of  absorbing  mid  scattering  material  and  the  same  tempera¬ 
ture  as  the  original  laser.  T  his  artificial  layering  has  been  found  necessary  ( Wollenweber. 
l9SSa)  to  remove  the  anomalous  dip  ( Hughes.  I9K7)  w  hich  occurs  when  aerosols  are 
included  in  the  I.OWTRAN  6  radiance  calculations  for  zenith  angles  close  to  90  degrees. 
Since  the  AG  A  scanner  could  not  be  accurately  plumbed,  the  zenith  angle  of  the  infrared 
horizon  (the  pixel  line  corresponding  to  the  maximum  radiance)  was  taken  to  be  one-half 
ot  a  pixel  line  less  than  the  angle  at  which  the  I.OWTRAN  calculations  indicated  the 
relracteit  ray  pat 1  first  struck  the  earth  (<•)  =  90.(79  degrees  on  29  September  and  90.174 
degrees  on  25  November). 

In  Fig.  3.  the  solid  lines  represent  the  loci  of  points  which  allow  the  I.OWTRAN 
calculations  to  match  the  measured  horizon  pixel  radiances  tor  different  combinations  of 
air  mass  lactor  and  v  isibility.  Phis  feature  of  the  calculations  results  from  the  visibility  scal¬ 
ing  (actor  of  the  size  distribution  remaining  nearlv  constant  for  any  appropriate  combina¬ 
tion  o!  the  two  factors  and  from  the  relative  insensitivity  of  the  calculated  extinction  coeffi¬ 
cients  for  the  tar  infrared  wavelengths  to  the  air  mass  factor.  At  the  time  of  the  IR  radiance 
measurements,  upwelling  solar  radiances  were  detected  by  Channel  1  (0.58  0.65  /um)  of  the 
WHRR  on  board  the  NOAA-9  satellite  as  it  passed  close  to  San  Diego.  These  radiances, 
used  to  determine  the  total  atmospheric  optical  depths  (  rn).  were  themselves  determined  by 
using  the  SMC  (Science  Applications  International  Corp.)  satellite  radiance  computer 
code,  which  makes  use  of  a  direct  linear  relationship  between  upwelling  radiance  and  the 
total  atmospheric  optical  depth  (Griggs.  1975).  Although  the  vertical  structure  of  the 
meteorological  parameters  which  control  aerosol  growth  were  different,  the  total  optical 
depths  lot  the  two  days  were  nearlv  identical.  In  a  manner  similar  to  the  calculations  of 
infrared  radiance,  calculations  ol  total  atmospheric  optical  depth  were  made  to  agree  with 
the  measurements  made  on  both  davs.  In  calculating  the  optical  depths,  it  was  assumed 
that  till  the  aerosols  were  confined  to  the  mixed  boundary  layer  as  determined  by  the 
aircraft  I  lights.  I  he  calculations  ol  optical  depths  with  different  combinations  of  air  mass 
lactor  and  visibility  are  also  shown  in  Fig.  3  ax  the  dashed  lines  and  are  seen  to  be 
extreme!}  sensitive  to  the  air  mass  factor.  I  he  intersections  of  the  solid  and  dashed  lines 
then  determine  the  best  combinations  of  air  mass  factor  and  visibility  to  be  used  in  the 
background  calculations  1  he  inset  in  the  figure  shows  the  good  agreement  between  the 
interred  visibilities  and  those  calculated  by  using  aerosol  size  distributions  measured  at  the 
lowest  levels  o|  the  aircraft  Rights.  I  hese  visibilities  seem  reasonable,  since  the  I  .os  Coro¬ 
nados  Islands,  located  between  25  and  30  km  off  the  coast  ol  San  Diego,  were  not  visible  to 
the  naked  eve  on  either  day. 
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Figure  3.  Loci  of  points  of  LOWTRAN  6  calculations  with  different 
combinations  of  air  mass  factors  and  visibilities  which  match  measured 
values  of  IR  horizon  radiances  (solid  lines)  and  satellite-detected  visible 
optical  depths  (dashed  lines)  for  29  September  and  25  November  1987. 

\\  h i I c  N  AM  was  initially  developed  by  using  surface  measurements  of  aerosol 
M/e  distributions  and  meteorological  parameters,  it  is  allowed  to  vary  with  altitude  in 
I  OVV  I  RA\  (i  according  to  relative  humidity  up  to  an  altitude  of  2  km,  at  which  point  the 
I  topospherie  Model  is  introduced.  A  dynamic  model  using  NAM  as  the  surface  kernel  (the 
N.iu  Ocean  Vertical  Aerosol  Model  (NOVAM)l  has  recently  been  introduced  (dc  l.eeuw 
ei  a!..  1989)  to  better  account  for  the  vertical  variations  in  si/e  distribution.  It  was  devel¬ 
oped  through  a  cooperative  effort  between  the  Naval  Research  Laboratory,  the  Naval 
Postgraduate  School,  and  the  Naval  Ocean  Systems  Center  (Gathman,  i989).  The  Physics 
and  I  lectromes  Laboratory  f  NO.  I  he  Netherlands,  recently  joined  in  the  evaluation  and 
retmement  el'lort  (dc  l.eeuw  cl  til..  1989).  The  model  is  based  on  the  physical  processes 
aliening  the  production,  mixing,  deposition,  and  si/e  of  the  aerosols  within  the  marine 
atmosphere.  Aerosols  with  similar  origins  are  represented  by  separate  log-normal  si/e 
distributions.  I  lie  net  optical  effect  produced  is  the  result  of  the  superposition  of  all  the 
,iei  osol  gr  oups.  I  o  determine  the  aerosol  si/e  distribution  at  any  particular  level,  one  of  a 
M-t  ol  mixing  profile  models  is  used  The  selection  process,  indicated  by  the  flow  chart  in 
I  lg  4.  is  determined  by  the  meteorological  input  data  available  and  the  wavelength  at 
which  calculations  are  to  be  made.  In  addition  to  the  general  profiles  of  temperature  and 
humiditv  as  determined  bv  radiosondes,  the  inputs  required  by  NOVAM  include  the 
nhlowing  sot  ol  surface  observations: 

•  24-b  uir  averaged  wind  speed  (m  s) 

•  \ii  mass  (actor  (  I  ...  10) 

•  (  loud  cover  (tenths) 

•  (  loud  tv  pe  (  I  .  .  .  10) 
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Figure  4.  Flow  chart  for  NOVAM. 


•  10. 6-/am  extinction  coefficient  ;!:m  ') 

•  Present  weather  in  standard  code  (0.  .  .99) 

•  Height  of  lowest  clouds  (m) 

•  Zonai  seasonal  category  ( I ..  .6) 

•  Sea  surface  temperature  (°C) 

•  Air  temperature  (°C) 

•  Relative  humidity  (f7) 

•  Optical  visibility  (km) 

•  l  ocal  wind  speed  (m  s) 


I  he  product  of  \OVAM  is  primarily  a  file  of  the  extinction  and  absorption  coeffi¬ 
cients  at  various  levels  in  the  marine  atmosphere.  While  NAM  has  undergone  numerous 
evaluations  and  considerable  improvement  through  the  use  of  transmission  and  sky 
radiance  measurements.  NOVAM  has  not  yet  received  a  critical  evaluation.  Such  evalua¬ 
tion  under  well-defined  environmental  conditions  is  needed  before  it  can  be  incorporated 
into  LOW  TRAN.  Data  for  preliminary  evaluations  of  NOVAM  were  obtained  from 
measurements  taken  on  the  upwind  side  of  San  Nicolas  Island  on  July  18.  1987,  and  have 
be  m  presented  elsewhere  by  de  l.eeuw  et  al.  ( 1989).  Figure  5  shows  examples  of  extinction 
profiles  at  both  visible  and  infrared  wavelengths  bands.  The  meteorological  input  data  w'ere 
obtained  from  a  tethered  balloon  platform  on  which  a  nephelometer  and  a  particle  measur¬ 
ing  system  (PMS)  particle  spectrometer  ( Knollenberg)  were  located.  Extinction  at  0.55  nm 
is  obtained  directly  from  the  nephelometer.  whereas  extinction  at  different  wavelengths 
may  be  calculated  from  the  si/e  distributions  obtained  from  the  PMS.  Figure  5a  is  the 
extinction  profile  for  0.55  ,um  measured  with  the  nephelometer  and  calculated  from  the 
measured  si/e  distributions  and  from  the  NOVAM  model.  The  data  show'  considerable 
scatter  in  the  measured  extinctions  with  altitude.  T  he  NOVAM  prediction  is  within  the 
envelope  of  the  scatter  better  than  750  of  the  time  for  this  pariicular  case. 

Figure  5b  shows  the  comparisons  between  the  extinction  coefficients  for  the  3  5-gm 
band  calculated  by  using  the  model  and  the  measured  si/e  distributions.  A  similar  compari¬ 
son  tor  the  8  I  2-/um  band  calculated  from  the  measured  si/e  distributions  and  the  NOVAM 
prediction  at  10.6  gmi  is  shown  in  F  ig.  5c.  The  present  version  of  NOVAM  underestimates 
the  extinction  in  the  infrared  bands  in  the  region  above  the  mixed  boundary  layer.  This 
may  be  a  result  ol  larger  particles  from  the  sea  surface  being  mixed  into  the  atmosphere 
above  the  apparent  temperature  inversion  through  the  process  of  entrainment. 

Additional  evaluations  of  NOVAM  have  also  been  undertaken  by  using  airborne 
data  acquired  in  the  Southern  California  area.  The  results  indicate  that  predicted  and 
measured  extinction  coefficient  profiles  have  similar  structure  but  may  differ  in  absolute 
value.  I  he  agreements  are  strongly  dependent  upon  the  proper  choices  of  air  mass  factors 
and  the  scaling  to  visibility  .  Default  profiles  of  humidity  and  extinction  do  not  in  general 
represent  the  measured  values. 
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EXTINCTION  AT  8-12  (1/km) 

(c) 

Figure  5.  Comparison  between  extinction  profiles  predicted  by  NOVAM  and  experimental 
extinction  profiles  derived  from  nephelometer  data  or  from  particle  size  distributions  using 
Mie  theory  (Knollenberg  data).  The  measurements  were  made  from  a  tethered  balloon. 

(a)  Comparison  of  the  NOVAM  prediction  for  a  wavelength  of  0.55  with  a  Mie  calculated 
profile  at  0.55  /im  and  the  extinction  profile  derived  from  the  nephelometer.  (b)  Comparison 
of  the  NOVAM  prediction  for  the  3-5-/tm  band,  (c)  Comparison  of  the  extinction  profile  as 
predicted  by  OVAM  for  10.6  /4m  with  the  Mie  calculated  average  extinction  profile  in  the 
8-l2-/im  band. 
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MARINE  INFRARED  BACKGROUNDS 


\ s  was  demonstrated  earlier,  the  LOW  I  RAN  codes  have  proven  to  he  a  versatile 
tool  in  predicting  atmospheric  radiance  above  the  hori/on.  They  do  not  yet.  however, 
include  the  capability  of  estimating  the  apparent  radiance  o]  the  sea  surface  for  different 
sea  states  and  meteorological  conditions.  A  model  has  recently  been  developed  of  the  effec¬ 
tive  radiance  of  the  sea  surface  as  a  function  of  the  viewing  zenith  angle  and  sensor  height 
i  W ollenvv eber,  |9X8b)  and  has  been  incorporated  into  LOWTRAN  6  tor  processing  on  an 
HP-9020  computer.  In  addition  to  the  path  emissions,  this  version  uses  the  Gaussian- 
distributed  wave  slope  model  of  Cox  and  Munk  (1954),  which  is  based  on  surface  wind 
speed  and  direction  relative  to  the  look  angle.  Cuncntly  iimited  to  NAM.  the  model  basi- 
callv  calculates  the  total  contributions  to  radiance  at  the  sensor  from  the  emissions  of  the 
intervening  atmosphere  and  the  sky  reflections  and  emissions  Irorn  the  wave  slope  surfaces. 

f  igures  6  and  7  show  the  comparisons  of  the  measured  and  calculated  infrared 
radiances  for  zenith  angles  within  about  1  degree  below  the  horizon  for  the  air  mass  factors 
and  visibilities  for  the  2  days  in  Tig.  3  (Hughes.  1 989a ).  In  both  cases,  the  major  contribu¬ 
tor  to  the  total  radiance  just  below  the  horizon  is  the  path  emission  .V( <*>)/, .  While  the 
ref  lected  sky  radiance  A  ( <-) )rvA  and  the  surface  emission  .V((-))v>  are  small,  their  contribution 
to  the  total  radiance  ,V(  (-))„„  at  this  low  level  (33  m)  of  observation  cannot  be  neglected.  It 
is  interesting  to  note  onlv  a  small  reversal  of  the  relative  magnitudes  of  the  reflected  sky 


Figure  6.  Comparison  of  the  measured  and 
calculated  IR  radiances  for  zenith  angles  about 
1  degree  below  the  horizon  by  using  the  air 
mass  factors  and  visibilities  determined  for  25 
November  1987. 
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Figure  7.  Comparison  of  the  measured  and 
calculated  IR  radiances  for  zenith  angles  about 
1  degree  below  the  horizon  by  using  the  air 
mass  factors  and  visibilities  determined  for  29 
September  1987. 


radiances  and  surface  emissions  between  the  two  sets  of  calculations,  which  demonstrates 
the  small  influence  of  the  wave  slopes  for  the  moderate  wind  speeds  on  29  September.  Both 
the  calculated  and  measured  total  radiances  on  both  days  are  in  good  agreement,  which 
places  confidence  in  the  model’s  usefulness  in  radiance  calculations  at  other  altitudes  and 
zenith  angles. 

The  selected  aerosol  models  were  used  to  calculate  the  contributions  of  the  path, 
sea,  and  reflected  sky  radiances  to  the  total  background  radiance  as  a  function  of  altitude 
and  zenith  angle.  In  Fig.  8  and  9.  examples  are  presented  of  calculations  for  both  days  with 
a  sensor  altitude  of  1000  meters.  For  zenith  angles  less  than  about  95  degrees,  the  major 
contribution  to  the  background  is  the  path  emission,  with  the  reflected  sky  radiance  being 
less  than  !0 ri  of  the  total  in  both  cases.  In  Fig.  10.  the  resulting  apparent  blackbody 
temperature  of  the  sea  versus  zenith  angle  is  compared  for  both  days.  Again,  the  higher 
apparent  temperature  for  zenith  angles  near  the  horizon  on  29  September  results  from  the 
path  emission  from  the  warmer  elevated  layers.  The  rapid  fall-off  of  path  emission  with 
increasing  zenith  angle  (i.c..  shorter  slant  paths  to  earth)  is  the  cause  of  the  decrease  in 
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Figure  8.  Calculations  of  contributions  of  the  path,  sea,  and 
reflected  sky  radiances  to  the  total  background  radiance  as  a 
function  of  zenith  angle  as  viewed  from  an  altitude  of  1000 
meters  on  25  November  1937. 


Figure  9.  Calculations  of  contributions  of  the  path,  sea,  and 
reflected  sky  radiances  to  the  total  background  radiance  as  a 
function  of  zenith  angle  as  viewed  from  an  altitude  of  1000 
meters  on  29  September  1987. 
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Figure  10.  Calculations  of  the  apparent  blackbody  tempera¬ 
ture  of  the  sea  versus  zenith  angle  as  viewed  from  an  altitude 
of  1000  meters  on  29  September  and  25  November  1987. 


apparent  temperature  on  this  day.  In  contrast,  the  temperature  increase  on  25  November 
is  the  result  of  the  increase  in  emission  from  the  sea  surface  (which  was  warmer  than  the 
air)  with  increasing  zenith  angles.  The  relative  contributions  of  the  three  components  will 
of  course  change  for  other  altitudes.  As  the  nadir  zenith  angle  is  approached,  the  appar¬ 
ent  temperatures  at  this  altitude  do  not  reach  the  measured  sea  surface  temperatures 
(  7\s  =  I5.2°C  on  29  September  and  =  17,6°C  on  25  November)  because  of  the  path 
emission  contributions  at  this  altitude. 

In  a  practical  sense,  the  sea  radiance  model  needs  to  be  improved  to  reduce 
computer  running  time.  In  its  present  form,  approximately  3  minutes  of  computational 
time  on  the  HP-9020  is  required  to  calculate  the  sea  radiance  at  one  zenith  angle  of  viewing 
as  compared  to  approximately  5  seconds  for  the  sky  radiance.  Approximately  36  minutes 
ol  computer  time  is  required  to  calculate  the  mean  sea  temperature  within  1  degree  of  the 
horizon.  The  major  contributor  to  the  computational  time  is  the  large  number  of  sky 
radiance  reflections  and  surface  emission  calculations  that  must  be  made  to  account  for  the 
radiance  from  all  the  wave  facets.  Computational  time  could  be  saved  if  the  surface  emis- 
sivity  and  reflectivity  calculations  were  made  at  a  single  representative  wavelength  instead 
ot  having  to  average  over  a  wavelength  band  for  each  wave  slope.  Another  way  to  reduce 
the  computational  time  is  to  develop  an  empirical  relationship  to  directly  relate  the  appar¬ 
ent  mean  sea  temperature  to  the  calculated  mean  sky  temperature. 
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To  investigate  this  possibility,  the  THF  RMOl  F.KNIX  system  was  used  to  deter¬ 
mine  the  mean  equivalent  blackbody  temperatures  eorresponding  to  an  area  1  degree  above 
and  1  degree  below  the  hori/on  in  thermograms  taken  during  different  surface  wind  speed 
conditions  (  Hughes.  1989b).  For  this  study.  18  thermograms  recorded  during  different 
meteorological  and  wind  speed  conditions  in  1987  and  1988  were  used  to  compare  the 
mean  sky  and  sea  temperatures.  In  12  of  the  eases,  the  actual  sea  surface  temperature 
measured  by  the  aircraft  radiation  thermometer  was  available.  For  the  full  data  set.  the 
mean  temperatures  were  well  correlated  (coefficient  of  correlation,  r-  0.92).  as  show  n  in 
I  ig.  II,  with  the  sea  tempeialuics  being  less  than  the  sky  temperatures,  as  indicated  by  the 
dashed  line  for  one-to-one  correspondence.  The  mean  temperatures  differed  the  most 
during  low  wind  speed  conditions.  These  differences  were  found  to  decrease  with  increasing 
wind.  Inherent  in  the  observations  are  the  actual  sea  surface  temperatures  and  wind  speeds. 
It  these  parameters  are  included  in  a  multiple  regression  analysis  for  three  independent 
variables,  the  following  relationship  is  obtained  for  the  mean  sea  surface  temperature. 

/'„,( sea),  in  terms  of  the  mean  sky  temperature.  r,„(skv).  the  current  wind  speed  (  F  ).  and 
the  actual  sea  surface  temperature  ( 

Tj sea)  =  1.09  7,„(sky)  +  0.37  T  +  0.24  /'„  10.5  (10) 

1  he  correlation  coefficient  for  the  multiple  regression  analysis  is  0.97.  In  this  case,  only 
12  samples  were  used  in  the  analysis  (as  compared  to  18  samples  in  the  linear  regression 
analysis  shown  in  f  ig.  I  1).  which  indicates  a  definite  improvement  in  the  correlation  when 
the  actual  sea  surface  temperatures  and  wind  speeds  are  included. 


Figure  11.  Regression  arratysis  comparison  of  the  measured 
mean  sky  temperatures,  Tm  (sky),  and  mean  sea  temperatures, 
f  m  (sea).  The  dashed  iine  corresponds  to  a  one-to-one  corre¬ 
spondence  between  the  two  parameters. 
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Earlier  work  (Hughes  and  Jensen,  1 988 )  has  demonstrated  that  the  sky  radiances 
measured  at  individual  pixel  lines  within  about  1  degree  of  the  horizon  can  be  modeled 
very  closely  by  using  the  modified  LOWTRAN  6  code.  Based  on  this,  the  replacement  of 
Tm(sky)  by  a  calculated  value  of  the  mean  sky  temperature,  T c( sky),  would  appear  justi¬ 
fied.  This  would  allow  the  mean  sea  temperature  to  be  directly  determined  from  Eq.  10 
and  measured  profiles  of  meteorological  parameters.  Additional  data  are  needed  to  verify 
that  such  a  relationship  for  predicting  the  mean  apparent  temperature  of  the  sea  is  mean¬ 
ingful  for  other  meteorological  conditions  sea  states,  and  altitudes  of  observation.  Such 
a  relationship  would  be  a  valuable  inclusion  into  a  system  designer's  handbook  as  a  tool 
for  predicting  infrared  sea  backgrounds  in  locales  where  only  meteorological  data  are 
available. 


AEROSOL  SENSING  WITH  LIDARS 

Measurements  of  slant-path  visibilities  are  required  as  inputs  to  the  LOWTRAN 
codes  for  scaling  the  selected  size  distributions  with  altitude.  Lidar  systems,  which  measure 
the  radiation  backscattered  to  a  receiver  by  aerosols  at  different  ranges  within  the  beam  of 
a  laser  pulse,  have  received  considerable  attention  in  recent  years  for  this  purpose.  The 
technique  of  inverting  a  single-ended  lidar  return  to  obtain  range-dependent  atmospheric 
extinction  coefficients  requires  an  assumption  concerning  the  relationship  between  the 
volumetric  backscattcr  and  extinction  coefficients.  Where  the  integration  is  performed  in 
the  forward  direction  (i.e.,  away  from  the  receiver),  the  extinction  coefficient  a(r0)  out  to 
the  range  where  the  transmitter  and  receiver  field-of-view  overlap,  rQ,  must  be  known.  If 
the  integration  is  performed  in  the  reverse  direction  (toward  the  receiver),  the  extinction 
coefficient  at  the  final  range,  o(r y),  must  be  known.  Ferguson  and  Stephens  (1983)  used  a 
novel  iterative  scheme  in  an  attempt  to  select  the  value  of  o(r^).  The  value  of  a(rf)  at  a 
close-in  range  (w  here  the  returned  signal  is  well  above  the  system  noise)  was  varied  until 
the  calculated  extinction  coefficients  as  a  function  of  range  a(r)  allowed  the  calculated  and 
measured  lidar  returns  with  range  to  agree.  The  chosen  value  of  o(rr)  was  then  used  as  a(r„) 
to  integrate  out  from  the  transmitter.  This  technique,  however,  is  limited  to  situations 
where  the  ratio  of  backscatter  to  extinction  is  a  known  constant  along  the  path. 

Solutions  to  the  single-scatter  lidar  equation  have  been  presented  for  reverse 
(Klett.  1985)  and  forward  (Bissonnette,  1986)  integration  where  the  relationship  between 
the  backscatter  and  extinction  coefficients  is  assumed  to  vary  with  range  according  to 

ft(r)  -  C(r)o(r)k  (II) 

where  k  is  a  constant.  For  forward  integration,  the  extinction  coefficient  as  a  function  of 
range  is  given  by 


air)  = 


exp  [Sir,,)] 
Ci r„)a(r„) 


C(r)]  exp  [,S(r)] 

2  J  [I  C(r)]  exp  [S(r')]dr' 


and  for  reverse  integration  by 


air) 


exp  [g(/y>j 

C(r  f)n(r  f) 


[I  C(r)]  exp  [5(r)J 

♦  2/r/[l  C(r)]  exp  [S(r')W 
r 


(12) 


(13) 
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where  the  constant  k  has  been  chosen  to  be  unity.  While  these  solutions  allow  for  variable 
backscatter  and  extinction  coefficients,  their  usefulness  requires  an  a  priori  knowledge  of 
the  ratio  C(r)  as  a  function  of  range.  Salemink  et  al.  ( 19X4)  determined  values  of  o  and  /? 
from  horizontal  lidar  shots  by  using  the  slope  method  when  the  atmosphere  appeared  to  be 
horizontally  homogeneous.  They  then  presented  a  parameterization  between  values  of  o 
and  relative  humidity  ( 3 3 ' \  §  RH  Is  87r; ).  When  the  parameterization  was  used  to  invert 
visible  wavelength  lidar  returns  in  the  vertical  direction,  the  derived  extinction  coefficient 
profiles  (using  radiosonde  measurements  of  relative  humidity  )  sometimes  agreed  reasonably 
well  with  those  measured  by  aircraft  mounted  extinction  meters.  In  contrast,  de  Leeuw  et 
al.  ( 19X6)  using  similar  ty  pes  of  lidar  measurements  did  not  observe  a  distinct  statistical 
relationship  between  backscatter  and  extinction  ratios  and  relative  humidity.  Fitzgerald 
(  19X4)  pointed  out  that  other  factors,  such  as  the  aerosol  properties,  can  strongly  affect  the 
relationship  between  (1  o  and  relative  humidity  and  that  the  power  law  relationship  of 
F.q.  I  I  is  not  necessarily  valid  for  relative  humidities  less  than  about  80rr-  A  unique  rela¬ 
tionship  between  C{r)  and  relative  humidity  which  is  dependent  on  the  air  mass  characteris¬ 
tics  is  vet  to  be  developed. 

A  lidar  inversion  algorithm  using  a  double-ended  lidar  technique  recently  developed 
independently  by  two  researchers  (Paulson.  1986:  Kunz.  1987).  In  this  technique,  the 
assumption  concerning  the  relationship  between  the  backscatter  and  extinction  coefficients 
is  eliminated  by  comparing  the  powers  returned  from  a  volume  common  to  each  of  the  two 
lidars  located  at  opposite  ends  of  the  propagation  path.  However,  the  receiver  gain  of  both 
lidars  must  be  accurately  know  n. 

The  mathematical  formulation  of  the  double-ended  technique  is  as  follows:  Consider 
two  lidars  separated  by  a  distance  d.  If  we  assign  the  origin  of  the  propagation  path  to  be 
located  at  lidar  1.  the  range-compensated  power.  S(r).  received  by  lidar  1  from  a  volume  at 
range  r  is  determined  by  the  single-scatter  lidar  equation  to  be 


Sir)  |  =  In  A',  +  In  [[Mr)] 


2  I  o(r')dr' 
'  o 


and  that  received  bv  lidar  2  is 


Sir)2  =  In  AS  +  In  [(Mr)] 


o{r')dr' 


(14) 


(15) 


Where  A,  and  AS  are  the  instrumentation  constants  for  each  of  the  lidars,  and  o(r)  and 
[Mr)  are  the  volumetric  extinction  and  backscatter  coefficients,  respectively.  If  the  scattering 
particles  are  assumed  to  be  spherical,  the  backscatter  coefficients  are  eliminated  by 
subtracting  Fq.  15  from  Eq.  14.  and  we  are  left  with 


o  '  r 


Since 


I  o{  r')dr  -  I  o[r') dr’  f  o[r')dr' 
r  o  o 


(17) 
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Equation  16  becomes 


S(r)l  S(r)-,  =  In  (  A,  A , )  4  J  a(r')dr'  +  2  /  *  a(r’)dr'  (18) 

<m  r 

Taking  the  derivative  of  Eq.  18  with  respect  to  range,  we  obtain 

o(r)  =  -  —  [A(r, )  5(r,)]  (19) 

4  dr 

The  determined  values  of  a(r)  can  then  be  used  with  the  system  constants  in  either  Eq.  14 
or  Eq.  15  to  determine  the  associated  backscatter  coefficients. 

The  double-ended  technique  has  been  used  recently  by  Hughes  and  Paulson  (1988) 
to  examine  the  effects  of  spatial  inhomogeneities  on  the  single-ended  lidar  inversion  algo¬ 
rithms.  Measurements  were  made  with  two  visioceilometer  (Lindberg  et  al.,  1984)  lidars 
over  a  0.9825-km  slant  path  near  the  ocean  on  the  Point  Loma  Peninsula  in  San  Diego, 
California.  Eidar  1  and  lidar  2  were  located  approximately  38  meters  and  135  meters  above 
mean  sea  level,  respectively.  Each  lidar  was  pointed  about  3  or  4  meters  away  from  the 
other's  receiver,  and  the  firings  were  offset  in  time  by  about  1  second  to  avoid  amplifier 
saturations  and  signal  contaminations.  The  data  presented  here  are  samples  of  two  returns 
taken  on  27  October  1986  during  a  period  of  reduced  visibility  and  when  conditions  along 
the  propagation  path  were  observed  to  be  varying  both  spatially  and  with  time. 

In  Fig.  12a.  the  S(r)  values  calculated  from  the  individual  lidar  returns  are  shown. 

In  each  case  the  ranges  are  referenced  to  the  location  of  lidar  1.  Similarities  in  the  larger 
atmospheric  irregularities  are  evident  in  each  data  set.  However,  there  are  differences  in  the 
fine  structure  of  the  individual  S(r)  curves  which  may  be  related  to  the  slightly  offset 
sampling  volumes  and  firing  times  of  the  two  lidars,  as  well  as  to  the  different  background 
scenes  viewed  by  each  lidar.  However,  the  background  only  adds  a  constant  value  to  the 
observed  backscatter  signals.  This  contribution  was  found  to  be  small  in  all  cases,  since  the 
signal-to-noise  ratio  for  both  lidars  exceeded  10  dB  at  a  range  of  800  meters.  In  any  event, 
data-smoothing  was  necessary,  and  an  1  1-point  (82.5-meter)  running  average  of  each  S(r) 
curve  was  determined  before  taking  the  differences.  The  differences  in  the  smoothed  S(r) 
curves  for  each  data  set  are  then  shown  in  Fig.  12b.  The  derivatives  of  the  S(r)  difference 
curves  were  calculated  with  a  running  range  interval  of  15-points  (1 12.5  meters),  and  the 
values  of  a(r).  calculated  from  Eq.  19,  are  shown  in  Fig.  13a. 

The  corresponding  backscatter  coefficients,  (Hr),  are  shown  in  Fig.  13  b.  They  were 
determined  by  using  the  calculated  values  of  a(r)  with  the  appropriate  5(r)  values  and  system 
constants  in  either  Eq.  14  or  15.  The  backscatter  coefficients  determined  for  lidar  1  and 
lidar  2  are  in  good  agreement.  The  slight  differences  occurring  within  a  range  of  450  meters 
may  reflect  the  precision  with  which  both  systems  could  be  calibrated.  Interestingly,  the 
backscatter  coefficients  do  not  show  the  striking  fluctuations  of  the  extinction  coefficients, 
which  leads  to  the  conclusion  that  C(r)  was  not  constant  along  the  path.  The  ratios  of  the 
calculated  extinction  coefficient  profiles  and  the  backscatter  coefficient  profiles  for  lidar  2, 
C(r).  were  determined  (assuming  the  value  of  k  in  Eq.  1 1  to  be  unity)  and  are  presented  in 
Fig.  13c. 
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M,  -  S(r) 


Figure  12.  Comparison  of  range-compensated  power  return,  S(r),  for 
(a)  lidar  1  and  2  and  (b)  their  differences  in  range  from  lidar  1. 
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(b) 


Figure  13.  Comparison  of  calculated  (a)  extinction  coefficient  a{r),  (b) 
backscatter  coefficient  0(r),  and  (c)  backscatter/extinction  ratio  C(r)  vs. 
range  from  lidar  1 . 
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I7  xt  i  net  ion  coefficients  were  calculated  as  ;■  function  of  range  ior  forward  and 
reverse  integration  by  using  the  individual  lidar  returns,  and  were  then  compared  with 
those  determined  bv  the  double-ended  technique  when  C(r)  was  assumed  to  be  either  a 
constant  or  allowed  to  vary  with  range  according  to  Fig.  13c.  I  he  results  ot  the  double- 
ended  lidar  technique  provide  the  close-in  or  far-away  boundary  values  tor  each  lidar. 

1  igurc  14a  shows  the  comparisons  lor  forward  integration  determined  when  C(r)  is 
assumed  to  be  constant  over  the  propagation  path.  There  is  little  or  no  agreement  between 
the  extinctions  determined  by  either  lidar  and  those  of  the  double-ended  measurements, 
even  though  the  boundarv  values  were  specified.  Near  a  range  of  600  meters,  w  here  the 
values  o!  Ctr  I  mcrease  (Fig.  13c).  the  extinctions  bv  lidar  I  exhibit  the  well-known  "insta¬ 
bilities''  m  that  they  tend  to  mcrease  without  bound  The  extinctions  by  lidar  2  also  exhibit 
the  instabilities  bv  tending  to  zero  at  the  farther  ranges.  These  data  demonstrate  the  sensi- 
tivitv  ol  the  instabilities  to  the  magnitudes  ol  the  boundarv  values. 

I  igure  14b  shows  the  comparisons  tor  reverse  integration  when  C(r)  is  assumed 
hi  be  constant.  Although  the  solutions  are  “stable."  there  is  little  agreement  between  the 
double-ended  measurements  and  the  extinctions  determined  bv  either  lidar. 

II  the  values  of  C(rt  are  allowed  to  vary  with  range,  as  shown  in  Fig.  13c,  excellent 
agreements  are  obtained  between  the  calculated  extinction  coefficients  of  lidars  1  and  2. 

I  he  extinction  coefficient  profiles  tor  both  the  forward  (Fig.  15a)  and  reverse  (Fig.  15b) 
integrations  are  nearly  identical  with  the  corresponding  double-ended  measurements  shown 
m  Fig.  I  3a. 

I  hese  results  demonstrate  that  if  the  value  of  Cir)  varies  with  range,  but  is  assumed 
to  be  a  constant,  neither  the  single-ended  forward  or  reverse  integration  algorithms  will 
allow  range-dependent  extinction  coefficients  to  be  determined  with  any  assured  degree  of 
accuracv  even  il  the  initial  boundarv  values  are  specified.  If.  however,  the  manner  in  which 
(in  vanes  is  specified,  both  the  forward  and  reverse  single-ended  inversions  for  this  data 
-et  reproduce  the  double-ended  measurements  remarkably  well.  Whether  the  same  is  true 
tor  other  situations  needs  to  be  determined. 

It  the  conditions  under  which  the  forward  inversion  algorithm  is  stable  can  be 
established,  then  a  single-ended  lidar  inversion  technique  would  be  possible  when  augmented 
with  a  close-m  measurement  of  extinction  and  measurements  to  relate  C(r)  to  air  mass 
characteristics  and  relative  humidity. 

While  the  works  ot  M ulders  <  1 984)  and  de  I  eeuw  et  al.  (1986)  have  concluded  that 
no  relationship  exists  between  Hr)  and  relative  humidity,  their  measurements  did  not 
account  tor  changes  in  the  air  mass  characteristics.  Whether  or  not  such  a  relationship  can 
evei  be  identified  in  a  practical  sense  is  yet  to  be  determined. 

In  situations  where  the  different  layers  ot  the  atmosphere  are  horizontally  homo¬ 
geneous.  i he  need  Im  knowing  the  relationship  between  the  backscatter  and  extinction 
coetlieients  can  be  eliminated  bv  comparing  the  range-compensated  powers  received  from 
eat  h  altitude  along  two  or  more  different  elevation  angles  ( Russel  and  I  iv  ingston.  1984; 
Paulson.  I9XU)  \ssummg  the  atmosphere  to  be  horizontally  homogeneous  with  extinction 
and  hack'catter  coetlieients.  o(/i)  and  fHli).  respectively .  which  vary  only  in  the  vertical 
dnection.  the  powers  receiver!  Irom  an  altitude  ,  along  a  slant  path  r,  elevated  <f>,  degrees 
and  along  a  slam  path  r ,  elevated  dq  degrees  (Fig.  16)  are  given  by 

6  /f(  h  ,  ) 

- exp 

(  /l  ;  S  III  (I1)  |  ) 


.  /i, 

2  I  n)h')(/h'  sin  d) 
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Figure  14.  Comparison  of  extinction  coefficients  a(r),  vs.  range  from 
lidar  1  calculated  from  lidars  1  and  2  by  using  (a)  forward  and  (b) 
reverse  integration  and  those  using  the  double-ended  technique  when 
C(r)  is  assumed  to  be  constant. 
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Pir.) 


Kfi(  h  | ) 
(/j,  sin  <£•>) 


-  exp 


,  A, 

2  /  o(h')dh'j  sin  (/>-, 
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where  K  is  the  system  constant.  Dividing  Eq.  20  by  Eq.  21,  we  get 


P(r,  )(h,  sin  <t>,)2 


exp 
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P(r2)(h |  sin  d>2)2 
or.  taking  the  logarithm  of  both  sides 
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Then 

~  •>* 

.  h 2  ,  h2  ,  h, 

t  -  I  a(h')dh'  =  /  o(h')dh'  i  o(h')dh' 

'/»,  0  0 

[,S(r,)  S(r4)]  [S(r,)  S(r,)] 

r„  =  -  (27) 

2(  I  sin  <£  I  sin  cf> , ) 

In  principle,  if  the  atmosphere  were  horizontally  homogeneous,  the  lidar  beam 
could  be  swept  in  elevation  and  an  incremented  profile  of  extinction  (and  consequently 
backscatter)  could  be  determined  by  using  the  returns  from  closely  separated  angles  (Kunz, 
1988).  The  smaller  angular  separations,  however,  place  stringent  requirements  on  the 
jirruracies  to  which  the  range-compensated  powers  must  be  measured.  If  we  assume  that 
the  percentage  errors.  5.  in  the  Sir)  values  are  equal,  the  percent  error,  A r„,  in  the  calcu¬ 
lated  optical  depths  is 

Ar,,  =  [(  r’(  -  r0)  t n  ]  »  100  (28) 
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where  r',  is  given  by 


{[.*>(/-,)(  1  ±6)]  [5(r4)(l±5]|  {[S(r , )( I ±<5)]  [.9(r-,)(  l±<5)]J 

T"  -  - : - (2S 

2(  1  sin  (t>2  I .  sin  0, ) 

Where  ( r(1)  is  a  maximum,  the  signs  of  6  are  such  that  the  products  of  S(r)(±5)  are  all 
of  the  same  sign;  then 


±[.S(r,)  +  S(/-4)  +  -S(r,)  +  S(r,)]6 

r  r  ,  -  - ; -  (30) 

2(  I  sin  </>-,  1  sin  d>, ) 

1  o  demonstrate  this  worst-case  condition,  the  extinction  coelficient  is  allowed  to  vary  line¬ 
arly  from  a  value  of  0.05  km  1  at  the  surface  to  1.0  km  1  at  a  height  of  1  km.  Assuming 
that  the  relationship  between  backscatter  and  extinction  were  invariant  with  altitude  (i.e., 
C(r)  -  C).  S(r)  values  were  calculated  for  heights  of  0.1  km  and  1  km  for  elevation  angles  of 
<t>  |  =60  degrees  and  =  30  degrees.  The  sensitivity  of  the  optical  depth  to  uncertainties  in 
the  .Sir)  valut,«  i,>  shown  in  Fig.  17  for  in(A'C)  =  0  and  5.4.  which  demonstrates  that  the 
errors  are  are  also  dependent  upon  the  relationship  between  backscatter  and  extinction.  In 
both  cases,  the  sensitivity  of  the  optical  depth  to  errors  in  S(r)  measurements  is  readily 
apparent. 


Figure  17.  Errors  in  optical  depths  as  a  function  of  errors  in 
S(r)  measurements  for  an  extinction  coefficient  profile  which 
varies  linearly  from  0.05  km-1  at  the  surface  to  1  km-1  at  a 
height  of  1  km  for  different  products  of  system  constant  K  and 
backscatter-to-extinction  ratio  C. 
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The  extent  to  which  the  atmosphere  is  horizontally  homogeneous  could  be 
determined  from  a  single  horizontal  lidar  shot  by  examining  the  gradient  of  the  range- 
compensated  power  return.  Then,  if  the  atmosphere  appeared  homogeneous,  i.e.,  a 
constant  negative  gradient  in  S(r).  the  optical  depth  between  any  two  altitudes  could  be 
detert.'ir?d  from  27.  tJowcvei.  tin.  within  the  ecnvectively  mixed  marine 

boundary  layer  rarely,  if  ever,  has  the  degree  of  homogeneity  required.  Paulson  ( 1989)  has 
recently  conducted  an  investigation  as  to  the  usefulness  of  the  two-angle  lidar  technique 
in  a  coastal  region.  In  these  studies,  data  were  taken  beneath  a  thin  stratus  cloud  layer  at 
about  500  meters.  Two  visioceilometer  lidars  were  operated  side-by-side  on  the  west  side 
of  the  Point  I.oma  Peninsula  and  pointed  west  overlooking  the  Pacific  Ocean.  A  series  of 
nearly  simultaneous  shots  were  made  with  this  arrangement.  Without  changing  the  orien¬ 
tation,  the  elevation  angle  of  lidar  1  was  increased  to  25  degrees  and  that  of  lidar  2  to 
50  degrees.  A  series  of  nearly  simultaneous  shots  were  also  made  with  this  configuration. 

Examples  of  the  horizontal  range-compensated  power  returns,  S(r),  for  a  5-point 
running  average  for  each  lidar  are  shown  in  Fig.  18.  The  two  lidar  returns  are  in  quite  good 
agreement,  but  the  irregularities  and  increasing  return  with  increasing  range  indicate  an 
inhomogeneous  condition.  The  S(r)  data  plotted  as  a  function  of  altitude  for  the  two-angle 
shots  are  show  n  in  Fig.  19.  The  effects  of  the  inhomogeneities  are  evidenced  by  the  fluctua¬ 
tions  in  the  S(r)  curve  for  lidar  1  above  and  below'  that  for  lidar  2  at  different  ranges.  The 
optical  depths  between  different  altitudes  determined  from  Eq.  27  are  shown  in  the  follow¬ 
ing  table: 


Table  1.  Optical  depths  calculated  from  different 
altitudes  up  to  a  maximum  altitude  of  475  meters 
on  17  May  1989. 


Lower  Altitude  (m) 

Optical  depth 

100 

0.81 1 

125 

0.437 

150 

0.584 

175 

0.597 

200 

0.647 

225 

0.584 

250 

0.688 

275 

0.150 

300 

0.260 

325 

0.260 

350 

0.342 

375 

0.492 

1  he  optical  depth  between  275  and  475  meters  is  only  0. 1 5,  while  that  from  375  to  475 
meters  is  more  that  three  times  greater  (0.49).  If  the  data  were  representative  of  a  horizon¬ 
tally  homogeneous  condition,  the  optical  depth  up  to  475  meters  should  consistently 
decrease  as  increases. 
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Figure  18.  Examples  of  measured  S(r)  values  obtained  from  parallel 
horizontal  lidar  shots  for  lidars  1  and  2  with  a  5-point  running  average. 


Figure  19.  Examples  of  two-elevation-angle  lidar  shots  (<^>2  =  25  degrees 
and  <t>2  =  50  degrees)  with  a  5-point  running  average. 


ASSESSMENT  OF  SYSTEM  PERFORMANCE 


The  computer  code  PREOS  (Performance  and  Range  for  EO  Systems)  presently 
resides  in  the  US  Navy’s  Tactical  Environmental  Support  System  (TESS).  The  present 
version  of  PREOS  addresses  the  ability  of  FL1R  systems  on  the  S3-A,  A-6,  and  P-3  opera¬ 
tional  aircraft  to  perform  range-dependent  tasks  such  as  the  detection,  classification, 
and  identification  of  surface  targets.  The  code  utilizes  polynomial  fits  calculated  with 
I.OWTRAN  3B  (Selby  et  al..  1976)  of  equivalent  gaseous  absorber  amounts  based  on 
meteorological  parameters.  The  aerosol  transmittance  model  in  PREOS  was  developed 
from  a  model  obtained  from  the  work  of  B.A.  Katz  at  the  Naval  Surface  Warfare  Center 
(see  Hughes  and  Richter.  1979).  This  model  employs  height-dependent  marine  and  conti¬ 
nental  components.  In  a  manner  similar  to  the  LOWTRAN  6  NAM,  the  marine  contribu¬ 
tion  is  scaled  to  surface  wind  speed,  and  both  contributions  contain  relative  humidity 
particle-size  growth  factors. 

The  agreement  between  surface  target  detection  ranges  as  reported  by  operational 
squadrons  and  those  predicted  by  the  code  has  been  disappointing.  However,  the  wide 
scatter  between  observed  and  predicted  ranges  (Fig.  20)  is  probably  more  characteristic  of 
operator  capability  and  the  quality  of  input  data  than  the  predictive  capability  of  the 
model.  The  meteorological  inputs  to  the  code  were  obtained  in  many  instances  from  radio¬ 
sonde  launches  which  may  have  occurred  several  hundred  kilometers  from  the  observation 
sites  and  may  have  been  separated  in  time  by  several  hours.  Predictions  of  performance 
range  are  also  critically  dependent  upon  the  actual  temperature  of  the  ship  and  of  its  natu¬ 
ral  background.  The  algorithms  in  PREOS  are  currently  based  on  a  fixed  temperature 
difference  between  a  rectangular  target  and  its  natural  background.  This  approach  neglects 
the  effects  of  a  wind-ruffled  sea  on  sky  reflections,  emissions  from  the  surface  wave  facets, 
and  contributions  of  the  intervening  atmosphere  to  the  total  background  radiance  scene, 
which  changes  with  viewing  angle  and  the  altitude  of  the  sensor.  Without  knowing  the 
ship's  temperature,  which  is  dependent  upon  its  history  (course,  speed,  and  surrounding 
meteorological  parameters),  it  is  questionable  that  the  accuracy  of  detection  range  against 
an  adversary  can  be  much  improved.  However,  a  ship  commander,  when  aware  of  his  own 
ship’s  past  and  future  courses,  can  use  the  prediction  algorithms  to  determine  the  ranges  at 
which  an  adversary  can  detect  or  track  the  ship  using  passive  1R  sensors.  These  stand-off 
ranges  are  of  primary  importance  in  estimating  the  time  allowable  for  evasive  actions 
against  guided  weaponry  launched  at  the  ship  or  for  the  deployment  of  countermeasures. 

Currently  available  computer  codes,  such  as  SIR  EOS  (Burns  et  al.,  1980)  (three- 
dimensional)  and  SIRS  (Batlev.  1978)  (two-dimensional),  are  capable  of  using  several 
hundred  individual  structural  elements  of  a  ship  to  model  its  composite  1R  signature.  These 
codes,  however,  are  quite  complex  and  require  extensive  running  times,  making  them 
impractical  for  shipboard  use  in  a  real-time  prediction  system.  A  modification  to  the  SIRS 
computer  code  has  recently  been  developed  at  the  Nava!  Surface  Weapons  Center.  This 
modification.  SHIPSIG  (Ostrowski  and  Wilson,  1985).  approximates  the  complex  struc¬ 
ture  of  a  ship  with  plane  elements  which  represent  the  ship's  temperature  at  zero-range  on 
an  average  basis.  For  a  given  v  iewing  direction,  the  simplest  representation  of  the  ship 
consists  of  a  single  vertical  and  horizontal  element,  with  the  observer's  orientation 
accounted  for  by  appropriate  area  components.  In  the  present  model,  the  horizontal  and 
vertical  elements  and  ship-stack  correction  factors  apply  to  a  guided  missile  frigate-class 
ship.  T  he  model  requires  as  inputs  the  ship’s  course  and  speed  as  a  function  of  time  from 
a  starting  geographic  latitude,  the  surface  wind  speed  and  direction,  visibility,  relative 
humidity,  air  temperature,  the  ship's  initial  temperature,  and  the  viewing  angle. 
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Figure  20.  Comparison  of  observed  FUR  detection  ranges  and  those  predicted 
by  PREOS. 


A  FI  IR  detection  algorithm  has  recently  been  developed  (Hughes,  1989c)  by  using 
the  SHIPSIG  model  and  the  sea  surface  background  model,  which  varies  with  sensor  alti¬ 
tude  and  viewing  (zenith)  angle.  As  neither  model  has  undergone  extensive  validation,  the 
algorithm  is  not  yet  incorporated  into  the  operational  version  of  PREOS.  To  demonstrate 
the  use  of  the  algorithm  as  a  tactical  decision  aid  to  predict  the  vulnerability  of  a  frigate- 
class  ship  to  detection  by  an  airborne  common  module  FLIR.  a  case  study  is  presented 
(using  the  actual  course  of  a  frigate  operating  off  the  coast  of  San  Diego.  California). 
During  a  5-hour  period  the  ship's  course  changed,  allowing  solar  heating  of  different  sides 
of  the  ship. 

f  or  this  study,  the  Piper  Navajo  aircraft  made  a  vertical  spiral  in  the  vicinity  of 
the  ship  to  obtain  the  profile  of  temperature,  relative  humidity,  and  pressure,  which  are 
required  inputs  to  the  modified  l.OWTR  AN  6  computer  code  for  calculating  the  sea  and 
skv  radiances.  A  Barnes  PRT-5  radiation  thermometer  was  also  onboard  the  aircraft  to 
measure  sea  surface  temperature  from  low-level  constant-altitude  flights;  the  temperature 
was  determined  to  be  I6.4°C.  The  vertical  profiles  of  temperature  and  relative  humidity, 
which  were  measured  at  1340  PS  I  approximately  9  km  oil  the  coast  ol  San  Diego.  Cali¬ 
fornia.  are  shown  in  Fig.  21. 
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Figure  21.  Profiles  of  air  temperature  and  relative  humidity  measured  with 
altitude  on  9  June  1988  off  the  coast  of  San  Diego,  California. 


Near  the  time  the  meteorological  parameters  were  obtained,  measurements  of 
8  12-/am  horizon  radiances  were  made  with  the  thermal  imaging  system  (AG A  THER¬ 
MOVISION,  model  780).  The  current  and  24-hour  averaged  wind  speeds  ( Vc  -  2.9  m/ s  and 
V  -  2.8  m  s)  and  the  vertical  profiles  of  meteorological  parameters  were  used  in  LOW- 
TRAN  6  to  calculate  the  radiance  that  matched  the  maximum  radiance  in  the  scene  for 
nonunique  combinations  of  air  mass  factors  and  visibilities.  At  the  time  of  the  measure¬ 
ments.  the  Los  Coronodos  coastal  islands  off  San  Diego  were  barely  visible  to  the  naked 
eye  at  ranges  between  25  and  35  km.  Choosing  a  visibility  of  37  km,  an  integer  value  of  3 
resulted  as  the  appropriate  air  mass  factor  for  the  LOWTRAN  6  NAM.  Figure  22  shows 
the  comparison  of  the  measured  and  calculated  IR  radiances  for  zenith  angles  within  about 
I  degree  above  and  1  degree  below  the  horizon,  with  an  air  mass  factor  of  3  and  a  visibility 
of  37  km.  Both  the  calculated  sky  (6  <90.17  degrees)  and  sea  (6  >90.17  degrees)  radiances 
are  in  good  agreement  w'ith  the  measured  values  for  this  low-wind-speed  case. 

The  selected  atmospheric  model  was  used  to  calculate  the  contributions  of  the 
propagation  path  and  sea  and  reflected  sky  radiances  to  the  total  background  radiance  as  a 
function  of  altitude  and  zenith  angle.  I  n  Fig.  23  the  total  apparent  blackbody  temperature 
ol  the  sea  background  from  the  three  contributors  is  plotted  versus  zenith  angle  for  sensor 
altitudes  of  500  meters  and  2000  meters.  At  the  500-meter  elevation,  the  dip  in  temperature 
at  about  97  degrees  is  a  result  of  the  rapid  fall-off  of  propagation  path  emission  with 
increasing  zenith  angle  (i.e..  shorter  slant  paths  to  earth  than  at  the  2000-meter  elevation). 
For  zenith  angles  greater  than  about  100  degrees,  there  is  little  difference  in  the  apparent 
temperatures  at  each  altitude,  and  both  approach  the  measured  sea  surface  temperature  at 
the  nadir  zenith  angle.  This  is  in  contrast  to  Fig.  10,  where  propagation  path  emissions 
between  the  sensor  and  sea  surface  made  significant  contributions  to  the  total  radiance  at 
the  nadir  zenith  angle  for  both  at  an  altitude  of  1000  meters. 

Figure  24  shows  the  course  of  the  guided  missile  frigate  USS  Brooke  ( FFGI)  off 
the  coast  of  San  Diego,  California,  on  9  June  1988;  this  course  was  chosen  to  demonstrate 
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Figure  22.  Comparison  of  measured  and  calculated  IR  radiances  for  zenith 
angles  above  (6  <  90.17  degrees)  and  below  (d  >  90.17  degrees)  the  horizon. 


Figure  23.  Total  apparent  biackbody  temperature  of  the  sea  background 
versus  zenith  angle  for  sensor  altitudes  of  500  meters  and  2000  meters. 
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Figure  24.  Course  of  the  USS  Brooke  (FFG1)  on  9  June  1988. 


the  model.  During  the  5-hour  time  period,  changes  in  the  ship’s  heading  allowed  solar  heat¬ 
ing  of  different  sides  of  the  ship.  As  the  ship  completed  the  course  and  returned  to  harbor, 
it  passed  close  to  the  AGA  thermal  imaging  system  located  near  the  entrance  to  the  harbor. 
Ihe  AGA  system's  data  processing  software  allows  subtraction  of  the  sea  background 
radiance  surrounding  the  ship  and  provides  a  histogram  of  the  temperature  distribution  of 
the  ship  pixels  within  the  chosen  rectangular  area.  The  mean  temperature  of  the  ship  (uncor¬ 
rected  for  atmospheric  effects)  was  19.7°C.  The  measured  radiance.  Y(meas),  of  the  ship  at 
a  range  r  is  related  to  its  actual  effective  blackbody  radiance,  /V(ship),  and  the  atmospheric 
emission,  Y(path).  along  the  path  by 

V(meas)  =  ,-V(ship)r(r)  +  .V(path)  (31) 

where  r(r)  is  the  atmospheric  transmittance  at  a  range  r.  The  range  to  the  ship,  determined 
by  using  the  known  vertical  dimensions  of  the  ship  and  their  angular  subtense  within  the 
field-of-view  of  the  AGA,  was  approximately  1.7  km.  L.OWTRAN  6  calculations  of  trans¬ 
mittance  and  path  emission  were  made  to  determine  the  temperature  equivalent  to  Y(ship). 

I  hese  calculations  resulted  in  an  adjusted  AGA  average  temperature  measurement  of 
20.5°C,  assuming  that  the  surface  emissivity  of  the  ship  was  unity. 

For  the  ship  model  calculations,  the  ship's  initial  position  was  taken  to  be  near  the 
entrance  to  San  Diego  harbor.  The  initial  ship  temperature,  the  ambient  temperature,  and 
the  relative  humidity  throughout  the  course  were  not  recorded  by  the  ship.  These  values 
were  taken  to  be  constants  as  measured  at  the  AGA  site.  The  surface  wind  was  south¬ 
westerly  (252  degrees  true),  and  the  depression  angle  of  viewing  was  essentially  broadside 
at  0.6  degree.  I  he  average  ship  temperatures  calculated  for  the  port  and  starboard  sides  of 
the  ship  as  a  function  ol  time  arc  shown  in  Fig.  25.  The  most  apparent  features  in  the 
temperature  responses  are  the  rapid  heating  of  the  port  side  and  the  gradual  cooling  of  the 
starboard  side  as  the  ship  steamed  westward  in  the  early  morning,  and  their  abrupt  cooling 
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Figure  25.  Average  temperature  of  port  and  starboard  sides  of 
the  USS  Brooke  (FFGl)  and  the  adjusted  AGA  measurements  as 
the  ship  entered  San  Diego  harbor. 


and  heating  after  1000  hours  following  the  southeasterly  course  change  at  0952  hours.  The 
magnitude  ol  the  port  side  average  temperature  is  approximately  1°C  greater  (for  an  emis- 
sivity  of  unity)  than  that  measured  by  the  AGA  system  as  the  ship  returned  to  harbor  near 
1245  hours.  II  indeed  the  emissivity.  e.  of  the  ship  was  0.9.  as  is  assumed  in  the  model,  the 
measured  average  temperature  would  be  in  better  agreement  (22.7°C).  Allowing  for  the 
uncertainties  in  the  meteorological  parameters  surrounding  the  ship  throughout  the  course, 
the  reasonable  agreement  between  the  adjusted  AGA  measurements  and  the  model  predic¬ 
tions  is  gratifying. 

In  the  performance  calculations,  the  l.OW  I  RAN  6  code  is  used  to  directly  calculate 
the  sum  of  the  ship  and  path  radiances  received  by  the  sensor  at  a  range  r  as 

V(r).v,  =  X(r  =  0),r(r)  +  .\\r)p  (32) 

where  V(r  -  0),  is  the  ship  radiance  at  zero  range  and  N{r)  is  the  path  radiance.  A'(r)A+/> 
is  then  converted  to  an  equivalent  blackbody  temperature.  T(r)^p.  by  an  iterative  solution 
to  Planck's  blackbody  formula.  Similarly,  an  equivalent  blackbody  temperature.  T(r)h.  of 
the  background  radiance  at  the  specified  altitude  and  range  is  calculated,  and  the  resulting 
apparent  temperature  difference.  XT(r)a  -  T(r)s.t,  Th(r),  is  determined.  The  range  at 
which  the  apparent  temperature  difference  is  equal  to  the  system's  minimum  detectable 
temperature  difference  ( M  DTD)  curve  determines  the  maximum  detectable  range  (MDR) 
of  the  ship.  I  he  IT.IR  system  MDTD  versus  range  (spatial  frequency)  curve  was  calculated 
hv  using  the  formulation  lor  a  hypothetical  ship  operating  against  a  rectangular  target.  In 
fig  26  and  27.  the  calculated  MDRs  for  the  USS  Brooke  by  an  airborne  FLIR  operating 
at  altitudes  ol  0.5  km  and  1.0  km.  respectively,  are  shown.  The  MDRs  could  have  been 
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Figure  26.  Calculated  MDR  envelopes  for  the  USS  Brooke  (FFGl)  by  an 
airborne  FUR  at  an  altitude  of  0.5  km. 
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Figure  27.  Calculated  MDR  envelopes  for  the  USS  Brooke  (FFGl)  by  an 
airborne  FUR  at  an  altitude  of  1.0  km. 


34 


calculated  as  a  hi  union  of  lime  throughout  the  the  ship's  course:  however,  lor  the  sake  of 
simplicity,  onh  the  v  ulnerabilitv  detection  envelopes  for  tire  entire  duration  ot  the  ship's 
course  are  shown  In  the  figures.  the  ship’s  aveiage  temperature  lor  both  the  port  and 
starboard  sides  are  shown.  In  Fig.  26.  the  ship  is  seen  to  he  v  ulnerable  to  !  R  detection 
throughout  its  course  Irotn  an  altitude  oi  0.5  km  at  a  range  of  31  km.  However,  beyond 
’3  km  the  ship  is  not  detectable.  Similarly,  in  Fig.  27.  the  ship  is  v  ulnerable  at  a  range  of 
s.'  km  from  an  altitude  ot  1.0  km.  but  is  safe  from  detection  beyond  56  km.  In  Fig.  2K  a 
comparison  is  shown  of  the  MDRs  calculated  with  the  current  algorithm  and  with  the 
method  that  assumes  a  constant  temperature  difference  between  the  ship  and  its  back¬ 
ground  ot  5  (.  When  the  current  method  is  used,  considerable  increases  (*&2()  km  at  an 
altitude  ot  2.0  km)  are  obtained  compared  to  the  fixed-temperature  method  in  predicted 
detection  ranees  with  altitude. 


DFTECTION  RANGE  (km) 

Figure  28.  Comparison  of  the  MDRs  (lower  envelope)  for  the  USS  Brooke 
(FFGi)  calculated  with  the  current  algorithm  and  with  the  method  that 
assumes  a  constant  temperature  difference  of  5°C  between  the  ship  and 
its  background. 


C  ONCLI  SIONS  AND  REC  OMMENDATIONS 


In  t h c  absence  of  transmission  or  radiometric  measurements,  the  LOW  I  RAN  codes 
must  currently  be  relied  on.  along  with  measured  meteorological  parameters  and  models  of 
aerosol  si/e  distributions,  to  predict  the  atmospheric  effects  on  EO  system  performance. 
These  codes  have  proven  to  be  a  versatile  tool  in  predicting  atmospheric  radiance.  NAM 
has  undergone  considerable  evaluation  and  should  be  considered  a  workable  model  in 
performance  predictions.  While  NOYAM  is  still  in  the  development  and  validation  stage,  it 
shows  promise  of  accounting  for  the  dynamic  changes  within  the  mixed  layer  and  should 
transition  into  I.OW  I  RAN  in  the  near  future.  Still  lacking  is  the  ability  to  measure  the  air 
mass  factors  and  slant-path  visibility  in  operational  conditions.  Emphasis  should  be  placed 
on  alternative  methods,  such  as  the  sky  radiance  technique,  rather  than  on  air  mass  charac¬ 
terization  for  determining  the  required  input  parameters. 

While  the  remote-sensing  potential  of  iidars  shows  promise,  the  capability  of  deriv¬ 
ing  extinction  coefficients  from  backscatter  measurements  has  not  been  demonstrated  with 
any  assured  degree  of  accuracy  when  the  atmosphere  is  inhomogeneous.  The  two-angle 
technique  certainly  would  be  useful  in  situations  where  the  atmosphere  can  be  shown  to  be 
horizontally  homogeneous,  provided  that  the  system  calibration  can  be  accurately  deter¬ 
mined.  I  he  double-ended  technique  resolves  the  requirement  for  knowing  the  relationship 
between  backscatter  and  extinction.  Elow’ever,  the  requirement  for  having  an  instrumented 
range  at  both  ends  is  impractical  in  an  operational  situation.  It  appears  that  the  double- 
ended  technique  is  best  suited  for  aerosol  studies  and  model  validation. 

The  reliability  of  the  sea  surface  radiance  model  that  uses  the  LOWTRAN  6  NAM 
to  accurately  represent  measured  values  for  low  wind  and  moderate  wind  speed  conditions 
has  been  demonstrated.  Whether  or  not  it  will  be  representative  of  other  wind  speed  condi¬ 
tions  needs  to  be  determined.  Also,  the  preliminary  evaluation  of  the  average  ship  tempera¬ 
ture  model  showed  promise,  since  this  model  responds  to  the  differing  solar  conditions, 
future  attempts  at  validation  must  ensure  the  accuracy  of  the  ambient  meteorolog.cal 
conditions.  Onboard  ground-truth  radiometry  measurements  of  the  temperatures  of  differ¬ 
ent  portions  of  the  ship  are  also  needed  to  aid  in  determining  the  accuracy  of  the  adjusted 
average  temperatures  inferred  from  the  AGA  measurements.  Finally,  a  controlled  experi¬ 
ment  with  an  airborne  operational  system  must  be  conducted  to  determine  the  validity  of 
the  predicted  detection  ranges  under  varying  meteorological  conditions. 
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